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Abstract: We report controlled manipulation of mode splitting in an 
optical microresonator coupled to two nanoprobes. It is demonstrated that, 
by controlling the positions of the nanoprobes, the split modes can be tuned 
simultaneously or individually and experience crossing or anti-crossing 
in frequency and linewidth. A tunable transition between standing wave 
mode and travelling wave mode is also observed. Underlying physics is 
discussed by developing a two-scatterer model which can be extended to 
multiple scatterers. Observed rich dynamics and tunability of split modes 
in a single microresonator will find immediate applications in optical 
sensing, opto-mechanics, filters and will provide a platform to study strong 
light-matter interactions in two-mode cavities. 
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1. Introduction 

Whispering-gallery-mode (WGM) optical microresonators with ultra-high quality factors and 
microscale mode volumes are of interest for a variety of scientific disciplines ranging from 
fundamental science to engineering physics. Significantly enhanced light-matter interactions 
HI make WGM resonators remarkably sensitive transducers for detecting perturbations in and 
around the resonator, e.g., virus/nanoparticle detection at single particle resolution J2l[3l|4) and 
ultrasensitive detection of micromechnical displacement J5J. Moreover, the coupling of opti- 
cal and mechanical modes mediated by enhanced radiation pressure within the microresonator 
provides a superb platform to study parametric oscillation instabilities and radiation pressure 
induced cooling of mechanical modes @. In addition, WGM microresonators with asymmetry 
(e.g., induced by structural deformations etc.) have been useful for investigating the correspon- 
dence between quasieigenstates and associated classical dynamics in mesoscopic systems (TJ. 
Level crossing have been demonstrated in microtoroids by tuning the microtoroid aspect ratio 
0. 

One interesting phenomenon associated with enhanced light-matter interactions is the split- 
ting of the initially degenerate cavity modes in the strong coupling regime [9 10l fTTl[T2ll . Mode 
splitting manifests itself as a doublet (two resonances) in the transmission spectrum of the res- 
onator. Here we show, for the first time, that by tuning the coupling strength using nanoprobes 
or subwavelength scatterers, the split modes can be manipulated individually or together form- 
ing a tunable two-mode microcavity. We demonstrate that the modes can be tuned to cross or 
anti-cross in frequency and linewidth. Particularly interesting experimental observation is the 
tunable transition between standing wave mode (SWM) to travelling wave mode (TWM) and 
vice versa. In order to explain the observed rich dynamics, we develop a model which takes into 
account the interactions between WGMs and finite number of sub-wavelength scatterers. Re- 
sults of the numerical simulations based on this model explains and confirms the experimental 
observations as well as suggests a rather surprising dynamics, i.e., the two split modes can be 
made to cross each other twice in linewidth with frequencies undergoing anticrossing in one of 
the linewidth crossing region. All the observed phenomena imply there are exceptional points 
in this cavity-scatterers system. 

2. Experiments, Theoretical Model and Discussions 

FigureQ]depicts the schematics of our experimental scheme which is composed of a WGM sil- 
ica microtoroid resonator coupled to two nanoprobes prepared by heat-and-pull of optical fibers 
on a hydrogen flame followed by buffered HF etching. To couple light in and out of the res- 
onator, a fiber taper is used. Positions of the nanoprobes and the fiber taper are finely controlled 
by 3D translational stages. We place the first nanoprobe in the resonator mode volume and fix 
its position when mode splitting is observed. We denote individual resonances of the doublet as 
ft)~(lower freuqency) and 0) + (higher frequency) modes with corresponding linewidths y~ and 
7 + (/ > 7 + X respectively. Then the second probe is introduced. This probe bends down and 
slides along the surface vertically as it contacts the rim of the microtoroid. Due to cone-like 
shape of the tip, vertical movement gradually increases its diameter allowing to simulate a scat- 
terer of increasing size within mode volume without changing lateral position. This does not 
cause any significant damage to the microtoroid as witnessed by no observable change in the 
value of Q factor. As the size increases, the nanoprobe starts disturbing the already established 



SWMs. The evolution of SWMs and the amount of disturbance applied to co^ and depend 
on the size and location of the second probe relative to the first. 



Fig. 1 . (a) Schematics of the experiment showing a fiber taper coupled to a microtoroid, and 
two nano-fiber tips introduced into the mode volume. Transmission spectra are captured by 
a photo detector (PD). (b) SEM image of a microtoroid. (c) SEM image of a fiber tip. Inset 
shows the enlarged image of the tip. (d) Cross-section of a microtoroid ring showing the 
position of a fiber tip in the field of a WGM. 

Figure [2] shows intensity graphs and transmission spectra recorded while the size of the 
second probe is increased. The second probe is used to tune the splitting between the modes 
and their linewidths. Crossing (Fig. |3b)), anti-crossing (Fig. |2jc)) and shifting (Fig. |2jd)) of 
linewidths and resonance frequencies at different lateral positions of the second probe are ob- 
served. 

In Fig. |2|b), the initial co + mode experiences red shift and linewidth broadening with in- 
creasing size of the second probe, while the co mode is not perturbed much. At a specific size 
both modes coincide, i.e., CO = co + and y~ = y + . Thus, a single resonance is seen in transmis- 
sion spectrum. With further increase of the second scatterer's size, the modes become separated 
with co mode now having a larger linewidth than co + whose linewidth equals to the initial y~. 
This suggests that both frequency and linewidth crossings have occurred. At the crossing point, 
back-scattering into the resonator vanishes as the back-scattered fields from the two scatterers 
have the same strength but ^-phase shift. Thus the vanishing of backward reflection coupled 
to the taper suggests a transition from SWM to TWM. The conditions for these to take place 
will become clear in the discussion of the theoretical model below. This observation implies 
that SWM, which limits nonclassical features of coherent matter-cavity field interaction due to 
the position dependence of the coupling strength, can be eliminated using external tuning with 
nanoprobes. This is particularly crucial in ultra-high-(2 microresonators because SWMs are 
usually formed due to mode splitting caused by structural defects and material inhomogenety. 

In Fig. Etc), the second probe first disturbs co + mode with no significant disturbance to 
CO . Thus co + experiences linewidth broadening and red shift gradually approaching to co~ . At 
a specific scatterer size, the frequency difference between the modes reduces from its initial 
value of 19.6 MHz to 7 MHz, and the linewidths become very close to each other. At this point, 
modes are strongly coupled to the scatterer and to each other. With further increase in size, CO 
strongly couples to the scatterer, and the modes start to repel each other leading to increased 
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Fig. 2. (a) Intensity graphs of mode- crossing (top), anti-crossing (middle) and shift (bot- 
tom). (b)-(d) Transmission spectra corresponding to the intensity graphs from top to bottom 
in (a). Increasing time corresponds to increasing size of the second nanotip. 



splitting. This suggests avoided-crossing of frequency and linewidth. 

In Fig. |2jd), the second probe affects both co and co + and induces frequency shift and 
linewidth broadening. The rate of change in CO is higher than that in co + suggesting that the 
scatterer has a greater overlap with co . 

The lateral position of the second scatterer determines whether the modes will undergo a 
frequency crossing, anti-crossing or shifting; however, the size of the second scatterer deter- 
mines the dynamics of the process. This is better understood by the following model: A single 
Rayleigh scatterer ifTTI with polarizability a,\ located at ri in the resonator mode volume V 
leads to a mode splitting quantified with the coupling strength 2g\ = —a\f 2 (r\)cOQ/V and the 
additional linewidth broadening 2T = a 2 f 2 (ri)cOQ /3nc 3 V where c is the speed of light, C0q is 
the resonance frequency before splitting, and / 2 (ri) is the spatial variation of the intensity of 
the initial WGM. The resulting two SWMs have periodic spatial distributions. A single scat- 
terer locates itself at the anti-node (node) of CO (co + ) with = (7t/2) where denotes the 
spatial phase difference between the first scatterer and the anti-node of an SWM El [Til . If a 
second Rayleigh scatterer with polarizability cti is introduced at location \2 with a spatial phase 
difference of j3 from the first scatterer, the already established SWMs redistribute themselves, 



and the amount of disturbance experienced by split modes depends on their overlap with the 
two scatterers (Fig.[3k). Subsequently, the frequency shift (Aco~ = ffl~ — (Oq) and the linewidth 
broadening (Ay~ = y~ — 70) of Q>~ with respect to the pre-scatterer resonance frequency coo 
and linewidth 70 become 

Acq- = 2gicos 2 (</>)+2 (?2 cos 2 (0-j3) (1) 
Ay" = 2ricos 2 0)+2r 2 cos 2 (> -]3) (2) 

where subscripts (1,2) represents the first and second scatterers. The cos 2 (-) terms scale the 
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Fig. 3. (a) Field distribution of a standing wave mode (SWM) obtained from finite-element 
simulation and the definitions of tj> and /3. Black circles represent the position of two scat- 
terers in the mode, (b) The relation between and /3. Solid curve is calculated from Eq. 
l[3} and * represent the values calculated using finite-element simulations for % = 0.5 and 
1 = 1. 

interaction strength depending on the position of the scatterer on SWMs [131 . Similar expres- 
sions for A© + and Ay + are obtained by replacing cos(-) with sin(-). Although we focus on 
two-scatterer case in this letter, this model can be extended to arbitrary number of scatter- 
ers by adding in Eqs. ([T]i and (f2| the terms 2g,-cos 2 (0 — J3/) and 2r,-cos 2 (0 — j3,), respectively, 
for each of the 3 < i < N scatterer. With the positions of the two scatterers fixed, the estab- 
lished orthonormal SWMs are distributed in such a way that the coupling rate between the two 
counter-propagating modes is maximized, in other words, the frequency splitting is maximized. 
It leads to: 

g 2 sin(2/3) sin(2j3) 



tan (20) 



gi+g 2 cos(2/3) £§ 2 + cos(2/J) 



(3) 



where ^ = f{r\)/f(t2) an d X = a i/ a 2 are positive real numbers. The two solutions of (j> have 
n/2 phase difference and correspond to the two orthogonal SWMs. Verification of Eq. <j3j is 
done by extensive finite-element simulations, and one example is presented in Fig.[3]5, where 
finite-element simulation results match very well with the calculated values from Eq. (0. 

Defining 5 = Aft) + — Aft) - = ft) + — ft) and p = Ay~ — Ay + as the frequency and linewidth 
differences of the resonance modes ft) and ft) + , we find 

5 = ^[Z 2 ^ 4 + 2^ 2 cos(2^) + l]3 (4) 

and 

p = ^^[x 3 <^ 4 +x(i+;^ 2 cos(2/3) + i] ( 5 ) 

from which frequency and linewidth crossings of the resonance modes can be calculated by 
setting 5 = and p = 0, respectively. 
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Fig. 4. Numerical simulations showing three unique patterns of doublet evolution for in- 
creasing size of the second nanotip at (a) /3 = n/2, (b)j8 = 0, (c)/3 = 0.44k for | = 1, (d) 
Resonance frequency and linewidth trajectories of the doublets with the exceptional point 
when % , x and are varied. Dashed and solid lines correspond to the two SWMs 



A. Behavior of the frequencies of the resonance modes: Conditions to observe frequency 
crossing is found by setting 8 = which implies J 2 ^ 4 + 2j^ 2 cos(2j8) + 1 = 0. It is satisfied 
only when cos (2/3) = — 1 or j3 = n/2. In the following discussions we only consider < /3 < 
n/2, as cos (2)3) is an even function and has period of n. Results of numerical simulations 
calculated from the model are shown in FiglH which coincide very well with experimental 
observations in Fig|2] Following are detailed discussions: 

(i) [} = 0. We find tan(2<j>) = 0, i.e., 0=0, implying that both particles locate at the anti- 
node of ft) - , and Ay~ is maximized. This leads to 8 = 2|gi|(l +X 2 )- Thus decreasing % 
increases 8 by pushing ft)~ further away from o + (Fig.|4fb)). Note that if the size of the second 
scatterer reaches above Rayleigh limit 03), it may start disturbing ft) + , too. 

(ii) j3 = n/2. The second particle stays at the anti-node of ft) + , thus increasing its size sig- 
nificantly affects the frequency and the linewidth of o + while its effect on co is minimal. 



Then we find 8 = 2|gi|(l — # -1 ^ -2 )- This implies only one frequency crossing which occurs 
at % = t, ~ 2 , i.e., g\ = g2- For £ = 1, frequency crossing occurs at % = 1 (Figs.@ta)). 

(iii) < /3 < tt/4. We have < cos2/3 < 1 (i.e., cosine is positive in the first quadrant 
of the unit circle) which implies that for a fixed j3 in this interval, 8 is always greater than 
zero (8 > 0) and it increases with decreasing %, that is with increasing a^. Physically, this is 
understood as follows. The second scatterer affects both SWMs with strengths depending on 
its overlap with each mode. In this interval of p\ the overlap of the second scatterer with co~ 
mode is always larger than that with co + . Consequently, as ctq. increases, co mode is further 
red-shifted increasing 8. 

(iv) 7r/4 < j3 <n/2. We have — 1 < cos2j3 < (i.e., cosine is negative in the second quadrant 
of the unit circle) implying that for a fixed j5 in this interval, 8 is always greater than zero 
(8 > 0); however, contrary to the case (iii) 8 has a minimum at x^ 2 = l/l cos (2j3)| > 1 with 
8mm = 2 \g\ 1(1 — X~ 2 ^~ 4 ) 1//2 > 0. The physical process is explained as follows. In this case, too, 
the second scatterer affects both SWMs with strengths depending on its overlap with each mode. 
When the size of the second scatterer is small, co + feels it strongly and undergoes frequency 
shift coming closer to CO as % decreases. This changes (j) and increases the overlap of the 
second scatterer with «~ leading to their stronger interaction which consequently, red-shifts 
CO and helps avoid crossing co + (Fig. [3c)). With a sufficiently large a,2, (j) > n/4- will be 
achieved which means CO will have larger overlap with the second scatterer than the first one. 

B. Behavior of the linewidths of the resonance modes: Setting p = in Eq. (0, we find the 
condition for linewidth crossing as l+<? 4 Z 3 + z(l+Z)^ 2 cos(2j3) =0 which can be satisfied 
only when cos(2j3) < or n/A < j3 < 7C/2, because both | and % are positive real numbers. 

(i) ft = n/2. We have = and the two scatterers locate themselves at the anti-nodes of the 
each SWMs, i.e., 1st scatterer at CO and 2nd scatterer at co + . Thus scatterers independently 
affect the two SWMs. We find that a linewidth crossing takes place at J 2 ^ 2 = 1, which means 
Y\ = F2- In this case if we also have % = § = 1, which gives g\ = g2 and F\ = Y2 implying the 
two SWMs have identical frequency and linewidth, but have orthogonal spatial distributions. 
The two SWMs merge to a TWM in the direction of the initial WGM. The other directional 
TWM vanishes as witnessed in experiments by vanishing of backward reflection in the fiber. 

(ii) 7T/4 < [5 < n/2. The roots of p = can be found by setting ^ 3 ^ 4 + %(l + %)^ 2 cos(2j3) + 
1=0. This is a transcendental third-order polynominal equation whose roots are too lengthy to 
give here. Given njA < /3 < n/2, either none or two positive real roots can be found for % at 
specific values of j3 and %. This suggests that two linewidth crossing points may be observed. 
Indeed the double crossing patterns are seen in the calculated patterns shown in Fig.|4ja,b). In 
both cases, one "symmetric" linewidth crossing (indicated by arrows in Fig.|4ja,b)) coincides 
with a frequency anti-crossing. From the plot of mode position 0, we see that the two SWMs 
"switch" distributions (<p shifts by 7t/2) around this point. This switching takes place around 
X^ 2 = 1 and is the source of the symmetry of linewidth crossing. The other linewidth crossing 
takes place around J 2 ^ 2 = 1, where the linewidth of one SWM changes significantly faster 
than that of the other one. This indicates that mode (red line in Fig. [J a) and blue line in Fig. 
@tb)) has much larger overlap with the second scatterer. Depending on whether £, > 1 or £, < 1, 
the "symmetric" linewidth crossing is observed before or after the other one. In experiments 
both scenarios were observed (Fig.|4jc,d)). 

Moreover, in the case that no positive real roots are found for p = 0, there is no linewidth 
crossing although one can always find % for specific ^ and j3 which minimize p. On either 
side of this minimum, p increases implying linewidth anti-crossing. This can be explained in a 
similar way as the frequency anti-crossing when 7r/4 < /3 < n/2. 

(iii) < /3 < n/A. In this case p ^ and similar to 8 in this regime, p increases as ct2 
increases. Neither crossing nor anti-crossing can be observed. 
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Fig. 5. Calculated frequency shift, linewidth broadening of the doublet and SWM position 
<(> as a function of the ratio of second scatterer size over the first one, for (a) j6 = 0.44;r, 
| = l/4, and (b) (5 = 0.447T, E, = 4. Dotted and solid lines correspond to the two SWMs. 
(c),(d) Experimental observations corresponding to (a) and (b), respectively. 



Note that the above discussed behaviors are associated with exceptional points (Fig. Hfd)) 
IT) and suggest that exceptional points can be observed in a single optical resonator by precise 
control of external nanoprobes. 

3. Conclusion 

In conclusion, we have discussed the controlled manipulation of WGMs in a microcavity using 
two external nanoprobes. The theory and method described here provide guidelines for manip- 
ulating the coupling of two SWMs and probing the mode splitting related phenomena. This 
work suggests that many application (e.g., filters, gyroscopes, delays) which are thought to re- 
quire multi-cavity configurations could be realized within a single microcavity. It also paves the 
ground for the discussions on tunable dual wavelength lasing with active material doped mi- 
croresonators, detection of consecutively introduced nanoparticles using microresonators and 
studying level dynamics of interacting modes. Continuous tunability of the split WGM modes 
by external probes will enable to selectively amplify mechanical modes as well as help efficient 
cooling of mechanical modes. 
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